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Note on dust trapping in inviscid vortex pairs with unequal strengths. 
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We investigate theoretically the motion of tiny heavy passive particles transported in a plane in- 
viscid flow consisting of two point vortices, in order to understand particle dispersion and trapping 
during vortex interaction. In spite of their large density, particles are not necessarily centrifugated 
away from vortices. It is observed that they can have various equilibrium positions in the reference 
frame rotating with the vortices, provided the particle response time and the vortex strength ratio lie 
in appropriate ranges. A stability analysis reveals that some of these points can be asymptotically 
stable, and can therefore trap particles released in their basin of attraction. A complete trapping dia- 
gram is derived, showing that any vortex pair can potentially become a dust trap, provided the vortex 
strength ratio is different from (single vortex) and -1 (translating symmetrical vortices). Trapping 
exists for both co-rotating or contra-rotating vortex pairs. In the latter case, particle trapping on a 
limit cycle is also observed, and confirmed by using Sapsis & Haller's method [Chaos, 20, 017515, 
2010] generalized to non-inertial reference frames. 
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I. INTRODUCTION 

Dust transport in vortical flows is a topic of interest in many research areas, since many natural 
or industrial fluid flows carry tiny solid objects. The understanding of particle dispersion, especially in 
turbulent flows, has been an active field of research in the last century (see for example Csanady^, Snyder 
& Lumley^i^, Maxey^). One of the most striking features of these particles, however, is their tendency to 
accumulate in various zones of the flow domain. This behavior has been observed also when particles are 
advected in a laminar flow : particles can be trapped by an attractor, whether chaotic or not, and remain 
there for long times (Rubin, Jones & Maxey^ ; Vilela & Motter~ ; Haller & Sapsis^). This makes the 
dynamics of inertial particles significantly different from the one of perfect tracers : in divergence-free 
flows the latter never accumulate. 



* Corresponding author : jean-regis.angilella@ensg.inpl-nancy.fr 



2 



Accumulation is rather unexpected for heavy particles in vortical flows. Indeed, these particles are 
generally centrifugated away from vortical zones under the effect of their inertia. In this paper we show 
that heavy particles can be trapped by vortex pairs, provided the particle response time and the vortex 
strength ratio lie in an appropriate range. 

Inviscid co-rotating point vortex pairs are known to rotate around each other under their mutual 
influence, and heavy particles released in such a flow do not have trivial trajectories. It was shown in 
a previous work (Angilellai) that particles released in a vortex pair with identical strengths can have 5 
equilibrium positions in the reference frame rotating with the vortex pair, and that two of these positions 
are asymptotically stable : they can attract particles. These equilibrium positions result from the balance 
between the hydrodynamic force, which drives the particle towards the center point of the vortex system, 
and the centrifugal force which thrusts the particles outward. In the following sections, particle motion 
in a vortex pair with unequal strength is examined for any vortex strength ratio. The flow is described in 
section im then particle trapping in equilibrium points is investigated (sectionllll]). In section|IV]we show 
that particles in counter-rotating vortices can also be trapped by a limit cycle. Particles with a very large 
response time are discussed in section FVl 




FIG. 1: Coordinate system and co-rotating frame (solid lines). The dashed axes are attached to the 
inertial reference frame. 

II. FLOW AND PARTICLE MOTION EQUATIONS IN THE CO-ROTATING FRAME 

We consider an incompressible and inviscid flow generated by a system of two point vortices with 
circulations Ti and located at Ai and A2. Elementary vortex dynamics shows that the distance 
2d between the vortices remains constant, and that they rotate with a constant angular velocity Q, = 
(Fi + /(Stt) around a point /, satisfying TiZZt + r2/^ = 0, provided Ti + r2 / 0. In the 
reference frame rotating with the vortex pair the flow is steady, and vortices remain fixed. In non- 
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dimensional form, using 1/0 for time and d for lengths, the vortex positions in this frame are : 



where the coordinate system (x, y) is centered at /, with Ix parallel to A1A2 (Fig. \Q. The fluid velocity 

dy 



field reads (x, 7) = ( §f ) — f§ ) 1 where ^' is the streamfunction in the co-rotating frame which can 



be written as : 

$(f) = - xi) + ^>(f - X2) + |x|^/2, 

with = — In |x| . Vorticity is constant and equals —2 everywhere, except at the vortex centers. 
Figs. [3] and m shows typical relative streamlines ^ = constant for various strength ratios 7. We will 
limit ourselves to 7 e] — 1, 1], since the case I7I > 1 can be treated by exchanging the role of vortices, 
i.e. by setting 7 — )• I/7. 

Particles are assumed to be non-interacting, non-brownian, much heavier than the fluid, and with a 
low Reynolds number (based on their radius and on the slip velocity). Clearly, the first approximation is 
no longer valid when particles tend to accumulate (see for example Medrano et al.-). However, we will 
use it for the sake of simplicity in this paper, and keep the effect of particle interactions as a perspective 
to this work. The simplest motion equation of these heavy particles at low Reynolds number reads, in 
non-dimensional form : 



5/ = '''(-^')-ir+n-^'-^'^''<^i <" 



where Xp{t) = {xp, Up) denotes the particle position at time t, Cz is the unit vector along the z axis, and r 
is the non-dimensional response-time of the particle (Stokes number). For spherical particles with mass 
nip and radius a the Stokes number reads r = Qnip / {<6t: fia) , where is the fluid viscosity. The last 
two terms in equation ([T]) are the (non-dimensional) centrifugal and Coriolis pseudo-forces respectively. 
Because particles aie much heavier than the fluid, Eq. ([T|) is valid even if the Stokes number r is not 
small, as added mass force, pressure gradient of the undisturbed flow, lift and Basset force are negligible. 
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III. TRAPPING CONDITION. 

Introducing a new vector variable z = (xp, y^, Xp, {jp), equation O can be rewritten as a system of 
first order equations: 

Zl = Z3 
Z2 = 24 

(2) 

Z3 = -Z3/T + Zi + 2Z4 + U^{zi,Z2,^)/t 
_ Zi = -Zi/r + Z2- 2Z3 + Uy{zl,Z2,-^)/T, 

where Vf = {ux,Uy). In vector form tiiis set of equations writes dz/dt = F{z,'y,T). To find the 
conditions of particle trapping, we search for attracting equilibrium points of (O : F{z, 7, r) = 0. If all 
the eigenvalues of VF have strictly negative real parts, then X^q is an asymptotically stable hyperbolic 
point. It can be shown (see AppendixjAjl that the eigenvalues of VF at the equilibrium point only depend 
on the jacobian J = Ux^xUy,y — Uy^xUx,y (the comma indicates spatial derivation) of the velocity field 
at this point, and the stability criterion is — < J < 1. For our velocity field the condition J < 1 is 
automatically satisfied. Let Xgq be an equihbrium point : 

f{Xeg,l, r) = Xeg + VfiXeg, 1)/t = 0. (3) 

By virtue of the implicit functions theorem, X^q depends smoothly on (7, r) unless det(V/ \ x^q ) = 0. 
Because det(V/ ) = J/t^ + 1, this condition coincides exactly with the boundary of the stability 
domain. Hence, the point X^q will exist and remain stable until (7, r) reaches some critical value defined 
by the system of equations: 

^ Xeq + Vf{Xeq,l)/r = 

JiXeq,7) + r^ = 0. 

This is a system of 3 equations for 4 scalar variables (Xeg,r, 7), which defines a family of curves in 
a four-dimensional space. The projection of these curves on the (r, 7)-plane splits it into subdomains, 
within which the number (and stability properties) of equilibrium points does not change. It is therefore 
sufficient to investigate the stability of equilibrium points for just one set of parameters (7, r) from each 
subdomain to get the full picture. 

To find these critical curves we took a grid covering 7 €] — 1, 1] and for each 7 solved the system 
dU) numerically. The resulting trapping diagram is shown in Figl2l for 7 g] — 1, 1] and < r < 0.6. 
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(Trapping also exists for lai^ger r : particles with a large response time ai^e discussed in the last section.) 
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FIG. 2: Particle trapping diagram. 



As expected, for 7 = 1 the critical value of r coincides with the one obtained analytically for iden- 
tical vortices^ : r = 2 — \/3. Numerical simulations of particle clouds have been run for the sets 
of parameters (A,B,C,D,E,F) maifed on figure |2l At t = particles were randomly distributed in a 
square [—3,3] x [—3,3]. Two types of initial velocities were considered : (i) no-slip initial condition 
{Xp{0) = Vf{Xp{0))), and (ii) particles at rest in the non-rotating frame {Xp{0) = —Q x Xp{0)). In 
these runs, we observed that for both (i) and (ii) the qualitative behavior of particles was the same, and 
the percentage of trapped particles did not change significantly. (Note however that particles with a large 
initial velocity might be located outside the 4-dimensional basin of attraction of the equilibrium point(s), 
and could therefore avoid trapping. ) Figs. |3] and |4] present the resulting particle clouds (black dots) at 
t = 30 in cases A, B, D, and at t = 500 for F (particles take much longer time to converge in this case), 
by using the no-slip initial condition. We observe that some particles converge to the attracting points 
predicted by the diagram, as expected. Grey dots, in Figs. ^2) and|H2), indicate the initial position of 
trapped particles (this is a cut of the basin of attraction). Clearly, the typical size of this basin is of order 
unity so that the percentage of trapped particles is non-negligible, about 15-20 % in cases A and D. We 
have checked that no trapping occurs in cases C and E, in agreement with the trapping diagram of Fig. |2l 

These results show that any vortex pair can trap particles, provided the vortex strength ratio 7 differs 
from - I or 0. 
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FIG. 3: Flow streamlines, particle clouds (black points) and equilibrium positions (white squares) for 
co-rotating vortices with different sets of parameters : A: 7 = 0.4, r = 0.1 ; B: 7 = 0.2, r = 0.2. Grey 
dots in Fig. (2) indicate the initial position of trapped particles in case A. 
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FIG. 4: Flow streamlines, particle clouds (black points) and equilibrium positions (white squares) for 
contra-rotating vortices with different sets of parameters : D: 7 = —0.1, t = 0.05 ; F: 7 = —0.8, r = 0.4. 
Figure F has been rescaled to account for the large radii of vortices trajectories when 7 is close to -1. 
Grey dots in Fig. (2) indicate the initial position of trapped particles in case D. 
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IV. PARTICLE TRAPPING INTO A LIMIT CYCLE. 

Figure m^l) suggests that, apart from the attracting points, a limit cycle could exist for small r and 
relatively small negative 7. The existence of a limit cycle can be predicted by using a criterion by S apsis 
and Haller-. To adapt this criterion to our equations we repeat the procedure from these authors, using an 
asymptotic expression of the particle velocity in powers of r. This procedure leads us to the following 
inertial equation : 

fix 

with Vi{Xp) = Xp — 2ez x Vj — VVf.Vf. With the use of this equation we get the following necessary 
condition for the existence of a limit cycle r-close to a flow streamline 5o (defined by y) = ^0) '■ 

I{So) = jvi-ndS = and j V{Vi ■ n) ■ ndS < 0, (5) 

So So 

where n is the outer normal to Sq. The former equation in (l5]l is a consequence of the existence of an 
invariant manifold for the particle dynamics r-close to Sq. The latter manifests the fact that this manifold 
is attractive. 

Now the problem is to find a value 'I'o such that /(5o(^'o)) = 0. This is a one-dimensional root- 
finding problem, which can be solved by applying the simplest bisection method, which does not require 
any additional information on /. The value of / at any given ^ is obtained by numerical integration along 
the corresponding streamline. For 7 = —0.1 this procedure converges to the value of ^'o ~ —0.761. The 
simulation caiTied out for particles with r = 0.01 confirms that the limit cycle indeed lies quite close to 
this streamline (Fig. |5]l. 



V. PARTICLES WITH A LARGE RESPONSE TIME 



Figure|6]shows the complete trapping diagram up to r = 10 : it suggests that particle trapping exists in 
the limit of large r. Indeed, the upper critical curves in this diagram have a vertical asymptote at 7 — )• 0, 
that is T 2/7^ (see Appendix |B]) : the grey zones are not bounded from above. Therefore, a stable 
equilibrium point can exist no matter how heavy the particles ai^e. Trapping is unexpected for particles 
with a large r because their trajectories are ballistic : they ai^e weekly influenced by the fluid. However, 
there is no contradiction here. The equilibrium point in coordinate- velocity space Zeq = {xeq, Ueq, 0, 0) 
is stable only locally and the basin of attraction for this point can be very small. Indeed, particles 



9 



1.5 
1.0 

0.5 

0.0 
-0.5 
-1.0 
-1.5 

-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 

X 

FIG. 5: Particle cloud at i = 50 for 7 = —0.1, r = 0.01 (grey points) compared with the predicted limit 
cycle (dashed black line). 

injected with a zero initial velocity (in the rotating frame) in the vicinity of the equilibrium point, will 
be slowly pushed to the attractor by the flow (with a rate of order l/r, as shown in Appendix iBt. In 
contrast, particles injected with a finite velocity will continue their inertial motion and will therefore 
avoid trapping. 

There are some points on the diagram where the derivative of the critical curves is no longer defined. 
These are the bifurcation values of 7, i.e. the values that turn the determinant of the system (|4| at 
the solution point to zero : det = 0, where g{p) is the right-hand side of the system (01, p = 

{xeq, Veq, t) IS an extended vector of parameters. Solving this equation for 7 together with the system 
dU), we get the following values : 71 ~ —0.50, 72 ^ 0.97. We can see from the diagram that contra- 
rotating vortices with 7 < 71 can (in principle) trap particles with any r, while for each 7 > 71 there 
exists a range of particles that are never trapped. 

The second point 72 corresponds to the appearance of two more branches of the critical curve at 
r ~ 3.56, but numerical simulations show that the percentage of trapped particles in these conditions is 
neglible (around 1%), so this part of the diagram is out of our interest. 
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VI. CONCLUSION 



We have shown that inviscid vortex pairs can permanently trap heavy particles, provided the vortex 
strength ratio and the particle response time lie in the grey zones of the trapping diagram shown in Fig. 
|2] Trapping can occur in the vicinity of various attracting points rotating with the vortex pair, where the 
centrifugal force balances the Stokes drag, for both co-rotating or contra-rotating vortices. The trapping 
diagram shows that any vortex pair can potentially become a dust trap, provided the vortex strength ratio 
is different from (single vortex) and -1 (translating symmetrical vortices). Contra-rotating vortices with 
7 < 7i — 0.5 can trap particles with any r. In contrast, for 7 > 71 there exists a range of particles that 
are never trapped. Numerical simulations show that the percentage of trapped particles with small or 
moderate r can reach tens of percents. For particles with a large inertia (r ^ 1) the convergence rate is 
very slow (~ 1/r) and the percentage of trapped particles is negligible. 

In the case of counter-rotating vortices with small r and relatively small 7, a limit cycle was found 
by using Sapsis & Haller's method^. A similar particle behavior has been observed in a cellular vortex 
flow model used for the analysis of the solar nebula (Tanga et alM-). 

An important issue is the effect of viscosity (which smoothes the fluid velocity field near the vortex 
cores and finally leads to vortex merging) on the trapping phenomenon. It could be natural to expect some 
sort of trapping provided the particle trapping time is much less than the characteristic time of vortex 
diffusion. Particle trapping in the viscous case has indeed been observed numerically^ for identical 
vortices. However, particle trapping by a viscous vortex pair with any strength ratio is still an open 
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question. 

Trapping leads to large collision rates and could be of interest to understand the detailed dynamics of 
various particle-laden flows where vortex pairing is known to occur. For example, anti-cyclonic vortices 
in protostellar disks are known to trap particles, with the help of the Coriolis force due to the disk rotation 
(Barge & Sommeria^, Tanga et al.—, Chavanis^). The analysis of the behavior of particles during the 
pairing of anticyclonic vortices could therefore be a topic of interest. 

Appendix A: Stability condition of tlie equilibrium points. 

Consider the characteristic equation of Q, taking into account the properties of the flow (vorticity 
equals -2, divergence equals 0): 

+ + (2 + l/r^) + 1 + J/r^ = 0, 

where A denotes any eigenvalue of VF. Using the change of variables ^ = A + l/2r, we obtain a 
biquadrate equation: 

16r^C^ + 8r2(4r2 - 1)C + ^t^{2t'^ + 2 J - 1) + 1 = 

which has the following roots: = (1 — 4r^)/4T^ it \/— J /r. Consider two possible cases: J < and 
J > 0. 

When J < 0, C± are both real and C+ > Now if C,\ < 0, all the roots C± are imaginary, 
so 5R[A] = — l/2r and the equilibrium point is stable. If at least > 0, then C+'^'' = ity''^ and 
max5R[A] = ^JOr ~ 1/2t. After some simple manipulations we find that maxK[A] > <J4> J < — r^. 

When J > the roots of the biquadrate equation are complex-conjugate: = 1/4t^ — libz\/j/r = 
a lb ib. Using trigonometric representation and the Moivre formula, we find that C± ^ = r^/^(cos <^/2 it 
i sin 99/2) and c!± = — C± \ where r = Va^ + &^ , cos if = a/r, so maxK[A] = \r^^'^ cosip/2\ — 1/2t. 
Using the fact that = (1 + cos(^)/2 we find that max3?[A] > <J4> J > 1. Combining 

these results, we obtain the stabihty criterion : the equilibrium point Xf.q is stable if and only if — < 

J(Xe„7) < 1. 
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Appendix B: Asymptotic behaviour for large t 

First we consider the case of finite 7 7^ 0,-1, taking e = lasa small parameter, and looking 

for the equiUbrium point in the following form: 

^eq = CxB"^, Veq = Cy£^. (Bl) 

Injecting this ansatz into the velocity field and keeping only the leading order terms in the limit e — )• 0, 
we obtain : 

^ Cys' (1 - (i±2ia±7!)) 
^j/ =^ 2 (7 - i) 

The equilibrium equations (|3]) then read (to leading order): 

c^^k+i ^1 + + = 0, 

2(7 - lh)e + CyE^ = 0. 

To ensure the consistency of these equations in the limit e — >■ we have to set m = 2, = 1. This leads 
to a system of two equations for Cx and Cy, which can be readily solved : 

Cy = -2(1 - 1/7), 

C, = 2(l-1/7)(1-B), 

with i? = (7 + 1)(7^ + l)/7^. Therefore, an equilibrium point of the form (|B1I) does exist for large 
T. To investigate its stability we calculate the Jacobian of the fluid velocity field at this point : J = 
(1 -^2) + 0(e). Using the expression for the roots of the characteristic equations from App. |A]and the 
Moivre formula, we obtain the following leading order estimates for the real part of the eigenvalues : 

\ -e/2, > 1. 

Evidently, max 3f?[A] < 0, so this equilibrium point is always (locally) stable for lai'ge r. The convergence 
rate is, however, very slow for these heavy particles and decays as t~^. 

The asymptotics above doesn't work for 7 — >■ 0. To investigate this case we consider 7 as a small 
parameter and search for the equilibrium point in the following form : x = Cx^"^, y = Cyj"^. The 
same method as above leads to e = Aj"- with 7i = 2 and m = 1, and the corresponding solution is 
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= 1± Vl - and Cy = 2 A. Clearly, C^, Cy are real only if yl < 1/2. The value A = 1/2 
corresponds to the upper critical curves in the (r, 7)-plane, that is r = 2/7^ (we have checked that the 
two upper curves of the numerical trapping diagram do coincide with this asymptote when 7 — > 0). When 
^ < 1/2 two equilibrium points exist, corresponding to C^; = 1 + \/l — AA'^ and Cx = 1 — Vl — 4^^. 
For the former point the roots of the biquadrate equation (Appendix E]) read : 

C± = -1 ± / 

1 + Vl - 4^2 

These roots are always real and negative, so K[A] = — This means that the point is locally stable, with 
a slow convergence rate (as in the case above). Similarly, one can check that the latter point is always 
unstable in the limit r — )• 00, with a divergence rate of order unity. 
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